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We measure the branching fractions of B° — > J/t^ry-'^ decays with the complete Belle data sample 
of 772 x 10 6 BB events collected at the T(45) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. The results for the branching fractions are: B(B° — > J/ipf]) = 
(12.3 ±i;| ±0.7) x 10~ 6 and B(B° -> J/ipn') < 7.4 x 10~ 6 at 90% confidence level. The 77-77' mixing 
angle is constrained to be less than 42.2° at 90% confidence level. 

PACS numbers: 13.20.Gd, 13.20.He, 13.25.Hw, 14.40. Nd, 14.40.Lb 



The neutral B meson decays B° — > J/ipr]('' are medi- 
ated by the b — > ccd transition as shown in Fig. [I] For 
such final states involving an 77 or 77' meson, it is conve- 
nient to consider flavor mixing of r] q and 77^ defined by 



Vq = —7={uu + dd), r)s 
v 2 



ss, 



(1) 



in analogy with the wave-functions of u and <f> for ideal 
mixing pQ. The wave-functions of the 77 and 77' are given 



b — * 




B° 



FIG. 1: Quark level diagram of the leading transition B° 
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If the ss component in Eq. ^ for B° — > J/ijjrf'} decays 
is negligible, the branching fractions are related to the 
77 — 77' mixing angle </> as: 

B(B° -> J/t/V 



B(B° -> J/V>r?) 



tan 



(3) 



Using the measured values of the 77 — 77' mixing angle 
0-40° and B(B° -> J/ipr)) = (9.5 ± 1.7 ± 0.8) x 
10~ 6 3 in Eq. the expected branching fraction for 
B° —> J/tpv' is about 6.7 x 10 -6 . Other predictions give 
B(B° -> J/tprf) = (4 - 8) x 1CT 6 and 77-77' mixing 
angle 37° < <j> < 50° 0H?]. Recently, B° s -> J/^ryW 
decays have been observed by Belle [5]. The branching 
fractions in Ref. [5] and this paper provide good inputs 
for model predictions ^[7j. The existing upper limit for 
the B° — > J jvbrj branching fraction is 6.3 x 1CT 5 [9]. 

In this paper, we report a measurement of B° — > J/ipr] 
and a search for B° —> J/yji]' decays [TO]. The results 
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are based on a data sample that contains 772 x 10 6 BB 
pairs, collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ (3.5 on 8 GeV) collider [IT] op- 
erating at the T(4S) resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters ( ACC) , a barrel-like ar- 
rangement of time-of-flight scintillation counters (TOF) , 
and an electromagnetic calorimeter (ECL) located inside 
a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail else- 
where [121 . 

The data sample used in this analysis was collected 
with two detector configurations. A 2.0 cm beampipe 
and a 3-layer silicon vertex detector were used for the first 
sample of 152 x 10 6 BB pairs, while a 1.5 cm beampipe, a 
4-layer silicon detector and a small-cell inner drift cham- 
ber were used to record the remaining 620 x 1Q 6 BB 
pairs rj3j. GEANT-based Monte-Carlo(MC) simulation 
program is used to model the response of the detector 
and to determine the reconstruction efficiencies. Simu- 
lated events are generated with the EvtGen program [13] . 
except for the B° — > ,// iprj' (rj' — > p°j) mode for which we 
use the QQ program [T3] since EvtGen does not properly 
generates the p° mass and angular distributions. 

Charged tracks are selected by using the impact pa- 
rameters relative to the interaction point: dr for the ra- 
dial direction and dz for the direction along the positron 
beam. Our requirements are dr < 1 cm and \dz\ < 5 
cm. Identification of e + and e~ from J/tp decay uses 
information from the ECL, the CDC (dE/dx) and the 
ACC. Identification of p + and candidates uses a track 
penetration depth and hit pattern in the KLM system. 
Charged pions are identified based on the information 
from the CDC (dE/dx), the TOF and the ACC. 

Photon candidates are selected from showers in the 
ECL, which are not associated with charged tracks, and 
an energy deposition of at least 50 MeV in the barrel re- 
gion or 100 MeV in the endcap region. A pair of photons 
with an invariant mass in the range 117.8 MeV/c 2 < 
M 77 < 150.2 MeV/c 2 is considered as a 7r candidate. 
This invariant mass region corresponds to a ± 3a inter- 
val around the ir° mass, where a is the mass resolution. 

We reconstruct J/ip mesons in the l + l~ decay chan- 
nel (I = e or fx). Any photon within 50 mrad of e + or 
e~ tracks is included as well. The invariant mass is re- 
quired to be within —0.15 GeV/c 2 < M ee ( 7 ) — mju < 
0.036 GeV/c 2 and -0.06 GeV/c 2 < M w - m J/v , < 
0.036 GeV/c 2 , where mj/^ denotes the nominal J/ip 
mass [16] . M ee ( 7 ) and are the reconstructed invari- 
ant mass of e + e _ (7) and p, + p~ , respectively. An asym- 
metric interval is used to include part of the radiative 
tails. 



r\ mesons are reconstructed in the 77 and tt + it 7r° 
final states. The mass ranges are 497 MeV/c 2 < 
M 77 < 590 MeV/c 2 and 520 MeV/c 2 < M„+„-„o < 
557 MeV/c 2 . In the 77 final state, candidates in which 
either of the daughter photons forms a tt° together with 
any other photon in the event are rejected. In the 77 
final state, we require | cos^| < 0.9, where 6 V is defined 
as the angle between the momentum of either of the pho- 
tons and the boost direction of the laboratory system in 
the rest frame of the r\. 

rj mesons are reconstructed in the 7/7r + 7r~ [r\ —¥ 77) 
and /9°7 final states. The mass ranges are 930 MeV/c 2 < 
M m+7T - < 967 MeV/c 2 , 920 MeV/c 2 < M p o 7 < 
980 MeV/c 2 , and 600 MeV/c 2 < M p o < 900 MeV/c 2 . 
In the /3°7 final state, we require cos (9,^ < 0.6 and 
|cos# p | < 0.8. Here 9^ is the angle between the pho- 
ton momentum and the opposite of the boost direction 
of the laboratory system in the rj rest frame. The other 
angle, 6 p , is the angle between the 7r + momentum and 
the boost direction of the laboratory system in the p rest 
frame. 

B mesons are reconstructed in the J/iprf'^ final state. 
Signal candidates are identified using two kinematic vari- 
ables defined in the T(4S) center-of-mass (CM) frame: 
the beam-energy constrained mass, M hc = -E^cam — Pb 
and the energy difference, AE = E^ — E\, eam , where pjj 
and -Eg are the momentum and energy of the B candi- 
date and -Ebeam is the run-dependent beam energy. To 
improve the momentum resolution, the masses of the se- 
lected 7T°, rf'} and J/ip candidates are constrained to 
their nominal masses using mass-constrained kinematic 
fits. In addition, vertex-constrained fits are applied to 
rj — > 7r + 7r~7T°, rj — > r)ir + TT~ and J/ip — > candidates. 
We retain events with M hc > 5.2 GeV/c 2 and \ AE\ < 0.2 
GeV. The signal peaks in the region defined by 5.27 
GeV/c 2 < M bc < 5.29 GeV/c 2 , \AE\ < 0.05 GeV (or 
-0.10 GeV < AE < 0.05 GeV for B° J/ipf](r] -> 77) 
only). 

For events with more than one B candidate, which are 
usually due to multiple rf'> candidates, the candidate 
with the minimum value from the mass- and vertex- 
constrained fit is chosen. 

The combinatorial background dominated by two-jet- 
like e + e~ — > qq(q = u,d,s,c) continuum is suppressed 
by requiring the ratio of second to zeroth Fox- Wolfram 
moments R 2 < 0.4 [17] , 

After the continuum suppression, the background is 
dominated by BB events with B — > J/ipX decays, where 
X denotes any final state. Using an MC sample of generic 
BB decays corresponding to 10 times larger than the 
data, with all known and expected B° — s- J/ipX de- 
cays, we study these backgrounds. The backgrounds 
from B° -» J/ipK and B° — > J/ipn cannot be highly 
suppressed and peak in the M& c distributions but not in 
AE. 

A figure-of-merit (FOM) method is used to opti- 
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mize the selection requirements. The maximal value of 
Ng/^/Ng + Nb is chosen for each variable, where iV s is the 
number of expected signal events and Nb is the number 
of background events. 

For each requirement, N s is estimated from signal MC 
simulation as 



B -> y 14* T|(y Y) 



B ^ yyiftCiciP) 



N s = N BS B{B° -> J/H)eB s 



(4) 



where N B g is the total number of BB pairs, 6.7 x 10~ 6 
is assumed for B(B° —> J/iprj'), e is the signal efficiency, 
and B sec is a product of the branching fractions for sec- 
ondary decays. 

The value of Nb is estimated from 



N b = R MC N b , 



(5) 



where Rmc is the fraction of B — > J /ipX MC events in 
the signal region, and N b is the number of data events in 
the AE sideband region. 

Signal yields and background levels are determined by 
fitting the AE distribution for candidates in the Mb c sig- 
nal region. For B° — > J/ipr], the AE distribution is fitted 
using a signal probability density function (PDF), which 
is a sum of a Crystal Ball function [T5] and two Gaussian 
functions. For B a —> J/ipr)', the signal PDF is a sum 
of three Gaussian functions. The background PDF is a 
second-order polynomial function. The signal shapes are 
determined from MC samples. The background levels are 
floated. 

We use the B+ -> J/ipK* + (K*+ -> K+tt q ) decay as 
a control sample to correct the difference between data 
and MC in the fitted mean and width of the AE signal 
peak. We require the helicity angle Ok* in the K* + — » 
K + tt° decay to be less than 90 degrees. Here Ok* is the 
angle between the 7r° momentum and the opposite of the 
boost direction of the laboratory system in the K* + rest 
frame. This requirement primarily selects events with a 
high momentum ir° and produces a control-sample AE 
distribution that is similar to that in our decay. 

The signal PDFs are modified based on the differences 
of mean and width between data and MC in the control 
sample. From a fit to the control sample, we find that 
the mean values are shifted by (—3.85 ± 0.13) MeV. The 
width in data is (1.11 ± 0.03) times wider than in MC 
simulation. 

We determine the signal yields by performing an un- 
binned extended maximum-likelihood fit to the candidate 
data events: 

e -(N s +N b ) *r 

C = Y[ [N s P s (AEi) + NbPbiAEi)], (6) 

i 

Here N is the total number of candidate events, 
P s (AEi) and P b (AEi) denote the signal and background 
AE PDFs, respectively, and i is the event index. 
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FIG. 2: AE distributions for the two decay modes, B° — > 
J/ipr/ and B° — > J/ipr/'. The plots on the top are B° — > J/ipr], 
in the left side - B° — > J/iprj(rj — > 77), in the right side - 
— > 7r + 7r 7r ). The plots on the bottom are 
B" -> J/i>T)', in the left side - B° -> J/ipv'W -> P°l), in 
the right side - B° — > J/ipr)'(r)' —¥ tity + tv~). The solid curves 
are results of the overall fit while the dashed curves are the 
background shapes. 



In the B° — > J/ifjr] mode, we fit the B° — > J /iptlill) 
and B° — > J/iprj(ir + ir~ir°) candidate samples simulta- 
neously with a common branching fraction. The fit 
gives the branching faction (12. Si^) x 10 -6 , which 
corresponds to signal yields of 77.9 and 29.8 for the 
77(77) and ?7(7r + 7r~7r ) modes, respectively. For the 
B° —> J/ipr)' mode, we also fit the B° — > J/tJnf(riir + ir~) 
and B° —> J/tpij' (p°j) candidate samples simultaneously 
with a common branching fraction. The fit gives the 
branching fraction (2.2±!'g) x 10~ 6 , which corresponds to 
signal yields of 5.5 and 5.2 for the r]'(p a -f) and ?/(?77r + 7r~) 
modes, respectively. The results of the fits to the data 
are shown in Fig. [2j 

A 1.4% systematic error comes from the uncertainty 
in the number of BB pairs. The systematic error due to 
tracking is 0.35% for each charged track. The systematic 
error from the pion identification requirement is deter- 
mined from a study of the D*+ -> D°ir+{D a -> K~-k + ) 
control sample. The systematic error from lepton identi- 
fication is obtained from a comparison between the data 
and MC for 77 — s- e + e _ //i + ^~ events. 

A 3.0% systematic error due to n° detection is deter- 
mined from a comparison of the data and MC ratios for a 
large sample of ?/ — > 7r + 7r _ 7r° and 77 — > 3ir° decays. Since 
V 77 is similar to ir° decay, we also assign a 3.0% 



5 



systematic error for r\ — ¥ 77 reconstruction. 

We use the previously mentioned control sample B + — » 
J/ipK*+(K*+ -> K+tt°) to obtain the systematic errors 
from the i?2 < 0.4 requirement and signal PDFs. A 
3.3% systematic error is obtained by comparing the data 
and MC. The systematic errors from signal PDFs are ob- 
tained by comparing the fit results for the cases when the 
fitting parameters are fixed from either MC or data and 
from changes that result from varying each parameter 
by one standard deviation. The errors on the branching 
fractions are taken from Ref. [T5] . 

The systematic errors are summarized in Tables |T] 
and [Tl] The total uncertainty is calculated by summing 
all individual uncertainties in quadrature. 

The statistical significances of the observed B° — > 
J /ipr] and B° — > J /ipr]' yields are 6.3a and 0.4cr, respec- 
tively. The significance is defined as ^/ —2hx(Lo / L max ) , 
where L max (L ) denotes the likelihood value at the max- 
imum (with the signal yield fixed at zero). The back- 
ground is floated and hence the systematic error is in- 
cluded in the significance calculation. 

An upper limit has been evaluated for the branching 
fraction of B° — > J/ipr/ at 90% confidence level because 
of the low significance. The Bayesian upper limits are 
obtained from: 

pN poo 

/ L(n)dn = 0.9 / L{n)dn, (7) 
Jo Jo 

where N denotes the signal yield. We use a modi- 
fied likelihood function to obtain a conservative upper 
limit, which is smeared with the systematic errors in 
the branching fraction. The upper limit for the branch- 
ing fraction of B° —> J /ipr]' at 90% confidence level is 
7.4 x 10~ 6 . 

TABLE I: Systematic uncertainties for B -> J /ipr] (%). The 
combined systematic error is 5.9%. 



77 source 


77 7T + 7T 7T° 


Number of BB events 


1.4 


1.4 


Tracking 


0.7 


1.4 


Lepton-ID 


2.6 


2.6 


Charged 7T-ID 




2.4 


r] — > 77, 7T° selection 


3.0 


3.0 


PDFs 


1.3 


1.3 


R 2 < 0.4 


3.3 


3.3 


B{J/ip -> e+e-,n + fj,-) 


1.0 


1.0 


B(i] -> 77) 


0.5 




B(T] -> TT+Tr-TT ) 




1.2 


Total 


5.7 


6.4 



The signal efficiencies and the branching fractions are 
listed in Table [ml 

In summary, we measure the branching fraction 
B{B° J/iprj) = (12.3 ±i;f ±0.7) x 10~ 6 . The first 



TABLE II: Systematic uncertainties for B — > J/ipT]' (%). 
The combined systematic error is 7.5%. 



7/ source 


777r + 7T 


P°7 


Number of BB events 


1.4 


1.4 


Tracking(lepton and charged pion) 


1.4 


1.4 


Lepton-ID 


2.6 


2.6 


Charged 7T-ID 


3.0 


3.0 


1 j.- 
r] — > 77, ty selection 


•3.0 


3.0 


PDFs 


3.6 


3.6 


R2 < 0.4 


3.3 


3.3 


B(J/ip -> e + e-,]i+]i-) 


1.0 


1.0 


B(r) -» 77) 


0.5 




B(v' -> A) 




2.0 


B(i]' — > r]n + n ) 


1.6 




Total 


7.5 


7.6 



TABLE III: MC efficiencies and branching fractions. 



Mode 




MC Efficicncy(%) 


/3(10- 6 ) 


B° — > 


J/iprjijj) 


35.7 




B°-> 


J/ipr](-K + ir~ 7T°) 


24.6 




B° -> 


J /ipr] combined 




12.3 ±h.f ±0.7 


B° -> 
B° — > 


J/iprj' (r]-K + ir~) 

J /Hip !) 

J /ipr] 1 combined 
J/ipr/ UL (90%) 


31.0 
19.1 


2.2 ±|;| ±0.2 
< 7.4 



error is statistical and the second is systematic. This re- 
sult is consistent with and supersedes the previous Belle 
measurement We do not observe a significant signal 
in B° — > J/iprj' and set the upper limit B(B° — > J /ipr]') 
< 7.4 x 10 -6 at 90% confidence level, which is eight times 
more stringent than the previous result [5] . From Eq. ^ 
we calculate the rj — r/' mixing angle, which is less than 
42.2° at 90% confidence level. These results are consis- 
tent with the theoretical predictions of Refs. [IHZ]- 
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